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ABSTRACT 

In  this  report  a  method  for  the  evaluation  of  cloud  temperature  and  transmittance 
from  ground-based  measurements  of  the  downwelling  spectral  radiance  is  proposed.  The 
method  uses  the  strong  emission  band  of  ozone  at  9.6  pm,  which  constitutes  a  natural  source 
of  IR  radiation  for  probing  clouds.  Temperature  and  transmittance  are  derived  from  the 
cloudy  sky  radiance  measured  in  two  narrow  spectral  channels  :  The  first  one  is 
approximately  at  the  center  of  the  ozone  band  at  9.48  pm  and  the  second  one  is  off  the  band 
at  9.15  pm.  The  cloud  parameters  are  found  by  solving  a  system  of  equations  requiring,  as 
input,  an  estimate  of  the  cloud  height  and  thickness  obtained  from  lidar  returns,  and  of 
atmospheric  temperature  and  humidity  profiles.  The  validation  tests  performed  on 
experimental  spectra  demonstrate  the  accuracy  of  the  method  with  typical  temperature  errors 
smaller  than  2  K  for  the  spectral  measurements  coincident  in  time  with  the  lidar 
measurements.  Finally,  the  main  sources  of  method  limitation  are  identified  and  analyzed. 


RESUME 

Dans  ce  rapport,  on  propose  une  mdthode  devaluation  de  la  temperature  et  de  la 
transmittance  des  nuages  en  se  basant  sur  des  mesures  de  radiance  spectrale  effectades  au 
niveau  du  sol.  Cette  mdthode  tire  avantage  du  fait  qu’il  y  a  une  forte  bande  d'dmission  de 
l'ozone  &  9.6  pm,  ce  qui  constitue  une  source  naturelle  de  rayonnement  infrarouge  pour  le 
sondage  des  nuages  La  temperature  et  la  transmittance  sont  dvaludes  k  partir  de  la  radiance 
du  ciel  nuageux  mesurde  dans  deux  regions  spectrales  dtroites:  la  premiere  se  trouve 
approximativement  au  centre  de  la  bande  d’ozone  k  9.48  pm  et  la  seconde,  k  l'extdrieur  de 
cette  bande  &  9.15  pm.  Les  paramdtres  du  nuage  sont  rdsolus  k  partir  d'dquations  qui 
ndcessitent  une  estimation  de  la  hauteur  et  de  l'dpaisseur  du  nuage,  et  des  profils  de 
tempdrature  et  d'humiditd.  Les  tests  de  validation  effectuds  sur  des  spectres  expdrimentaux 
ddmontrent  la  precision  de  la  mdthode  proposde,  avec  des  erreurs  typiques  de  tempdrature 
plus  petites  que  2  K  dans  le  cas  ou  les  mesures  spectrales  sont  en  coincidence  temporelle  avec 
les  mesures  du  lidar.  Finalement,  les  principales  sources  de  limitation  de  la  mdthode  sont 
identifides  et  analysdes. 
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EXECUTIVE  SUMMARY 


The  Canadian  Forces  are  using  an  increasing  number  of  electro-optical  systems  for 
the  remote  sensing  of  target  characteristics  and  other  information  of  military  interest. 
Defence  Research  Establishment  Valcartier  (DREV)  has  been  involved  in  remote  sensing 
for  many  years,  especially  in  the  field  of  lidars  and  Fourier  Transform  Infra-Red  (FTIR) 
spectrometers.  Recently,  DREV  has  demonstrated  a  new  multiple-field-of-view  lidar 
technique  for  the  remote  sensing  of  atmospheric  aerosols.  DREV  is  also  developing  a  new 
FTIR  technique  for  the  remote  sensing  of  atmospheric  parameters,  such  as  temperature  and 
humidity  profiles. 

One  application  for  which  these  two  DREV  techniques  appear  particularly 
promising  is  the  remote  sensing  of  icing  hazard  conditions.  The  military  and  commercial 
aviation  communities  are  very  concerned  with  this  problem  since  more  than  thirty-five 
aircraft  accidents  are  reported  each  year  in  North  America.  Aircraft-icing  conditions  happen 
in  a  cloud  composed  of  large  water  droplets  at  temperature  below  freezing  point.  In  this 
context,  we  have  recently  started  a  feasibility  study  of  an  on-board  lidar-radiometer  system 
for  mapping  icing  conditions.  The  first  step  of  this  project  consisted  in  performing 
simultaneous  lidar  and  spectral  measurements  (FTIR)  to  develop  a  method  for  determining 
cloud  temperature  and  composition:  ice  crystals  or  supercooled  water  droplets.  Lidar  data 
can  provide  information  on  particle  size  and  concentration  while  spectral  data  contains 
information  on  cloud  temperature. 

The  purpose  of  this  report  is  to  present  a  new  method  for  the  remote  determination 
of  cloud  temperature  and  transmittance.  The  method  is  founded  on  passive  measurements 
of  the  downwelling  spectral  radiance  coming  from  the  cloud.  This  development  is 
supported  by  ground-based  observations  performed  with  a  FTIR  spectrometer  having  a 
spectral  resolution  of  1  cm-1  and  covering  the  infrared  (IR)  band  from  3  to  20  (xm.  The 
cloud  evaluation  method  uses  the  strong  emission  band  of  ozone  in  the  IR  spectral  range. 
Because  most  of  the  ozone  emission  originates  from  the  upper  atmosphere,  it  constitutes  a 
natural  source  of  IR  radiation  for  probing  clouds  in  the  lower  part  of  the  atmosphere  using  a 
ground-based  receiver.  The  solution  for  cloud  temperature  and  transmittance  is  found  by 
solving  a  system  of  equations  that  requires  an  estimate  of  cloud  height  and  thickness  obtained 
from  lidar  returns,  and  an  estimate  of  atmospheric  temperature  and  humidity  profiles  for 
calculations  of  the  atmospheric  background.  The  method  has  been  successfully  tested  on 
experimental  spectra.  It  was  found  that  temperature  errors  are  typically  smaller  than  2°  C 
when  spectral  measurements  are  coincident  in  time  with  the  lidar  measurements.  The  main 
sources  of  errors  have  been  identified  and  their  impact  on  the  accuracy  of  the  method  has 
been  addressed.  The  next  steps  will  focus  on  generalizing  the  cloud  evaluation  method  to 
any  slant  paths  from  the  ground  or  from  an  airborne  platform. 

This  report  presents  a  basis  to  further  develop  the  cloud  evaluation  method  for 
applications  in  icing  hazard  detection.  The  ultimate  goal  is  to  demonstrate  that  combined 
lidar  and  spectral  measurements  (two  bands)  can  be  used  for  determining  the  temperature  of 
a  cloud  and  whether  it  is  made  of  ice  crystals  or  super-cooled  water  droplets.  The  results  of 
this  work  may  also  apply  to  other  domains  in  which  remote  sensing  techniques  are  utilized 
such  as  in  cloud  microphysics,  climatology  and  meteorology.  Finally,  this  work 
contributes  directly  to  Canadian  inputs  to  TTCP  JTP-13. 
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1.0  INTRODUCTION 

The  Canadian  Forces  are  using  an  increasing  number  of  passive  and  active  electro- 
optical  (EO)  sensors  to  assist  military  units  in  performing  tasks  such  as  surveillance, 
detection,  identification  and  designation.  The  continuous  effort  for  improving  sensors' 
performance  has  favored  the  development  of  many  sectors  of  science  and  technology.  One 
sector  that  has  experienced  very  rapid  and  successful  developments  is  remote  sensing.  For 
many  years.  Defence  Research  Establishment  Valcartier  (DREV)  has  been  deeply  involved 
in  this  field  of  research,  especially  in  the  development  and  exploitation  of  two  specific 
technologies:  lidar  and  Fourier  Transform  Infra-Red  (FTIR)  spectrometer.  In  the  case  of 
lidar,  DREV  has  recently  proposed  and  demonstrated  a  new  multiple-field-of-view 
(MFOV)  lidar  technique  for  the  remote  sensing  of  atmospheric  aerosols  (Ref.  1).  In  the 
case  of  an  FTIR  spectrometer,  DREV  has  been  developing  a  new  technique  for  the  remote 
sensing  of  atmospheric  parameters  such  as  temperature  and  humidity  profiles  (Ref.  2). 

One  application  for  which  these  two  techniques  (MFOV  lidar  and  FTIR 
spectrometer)  appear  particularly  promising  is  for  the  remote  sensing  of  icing  hazard 
conditions  for  aircrafts.  The  military  and  commercial  aviation  communities  are  very 
concerned  with  this  problem:  In  North  America,  more  than  400  aircraft  accidents  due  to 
icing  have  been  reported  since  1975  (Ref.  3).  Basically,  aircraft-icing  conditions  happen  in 
a  cloud  composed  of  large  water  droplets  found  at  temperature  below  freezing  point.  In 
this  context,  DREV  has  initiated  a  feasibility  study  of  an  on-board  lidar-radiometer  system 
for  mapping  icing  conditions.  The  first  step  of  this  project  consists  in  performing  a  series 
of  simultaneous  lidar  and  spectral  measurements  (FTIR)  to  develop  a  method  for 
determining  cloud  temperature  and  composition,  i.e.  ice  crystals  or  supercooled  water 
droplets.  Lidar  data  can  provide  information  on  the  particles’  size  and  composition,  while 
spectral  data  contain  information  on  cloud  temperature. 

The  purpose  of  this  report  is  to  develop  and  verify  a  new  method  for  the  remote 
determination  of  cloud  temperature.  The  method  is  founded  on  ground-based 
measurements  of  the  downwelling  spectral  radiance  coming  from  the  cloud.  In  Chapter  2 
the  FTIR  spectrometer  is  described,  while  Chapter  3  presents  recent  field  measurements  of 
the  downwelling  spectral  radiance  for  different  cloud  conditions.  It  is  shown  that  the  FTIR 
technique  is  particularly  efficient  in  observing  cloud  emission  and  scattering  because  of  its 
high  sensitivity,  very  wide  spectral  coverage  and  spectral  resolution.  In  Chapter  4  the 
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cloud  evaluation  method  is  formulated:  The  cloud  temperature  and  IR  transmittance  (9.3 
|im)  are  derived  from  the  cloudy  sky  radiance  measured  in  two  narrow  spectral  channels. 
Theoretical  and  experimental  verifications  of  the  method  are  reported  in  Chapter  5  together 
with  an  evaluation  of  the  two  main  sources  of  errors.  Finally,  Chapter  6  summarizes  the 
work  and  draws  conclusions  on  the  perspective  of  extending  the  cloud  evaluation  method 
to  scenarios  involving  airborne  measurements. 

This  work  was  performed  at  DREV  between  June  1995  and  May  1996  under 
Thrust  3.D  -  Airborne  Search  and  Localization,  Work  Unit  03D03-15:  Icing  Hazard 
Detection. 


2.0  DESCRIPTION  OF  THE  SPECTROMETER  SYSTEM. 

The  method  for  remotely  evaluating  cloud  temperature  and  transmittance  is  founded 
on  the  ground-based  measurement  of  the  IR  radiation  emitted  by  the  cloud  and  the 
surrounding  atmosphere.  For  the  present  study,  the  ground-based  measurements  were  taken 
with  a  DREV-designed  Fourier  spectrometer  referred  to  as  the  Double  Beam  Interferometer 
Sounder  (DBIS). 

Essentially,  the  DBIS  is  made  of  one  or,  optionally,  two  10-in.  diameter  Cassegrain 
telescopes  optically  coupled  to  a  double-input  port  Fourier  transform  spectrometer  and  two 
detection  units  (output  optics  1  and  2).  Figure  1  summarizes  the  design  of  the  instrument. 
Note  that  only  one  telescope  assembly  is  explicitly  shown  in  Fig.  1.  This  instrument  allows 
for  measurements  of  calibrated  spectra  according  to  the  following  specifications:  any 
selectable  zenith  angle,  scene  FOV  of  5  mrad,  spectral  coverage  from  3  to  20  Jim  and  a 
spectral  resolution  of  1  cm- 1  or  greater.  As  part  of  the  input  modules,  a  large  flat  plate 
mirror  placed  in  front  of  the  telescope  can  be  either  rotated  to  the  selected  scene  or  oriented 
in  a  position  for  the  acquisition  of  reference  spectra.  The  pointing  capability  of  this  scene 
mirror  makes  it  possible  to  take  slant  path  measurements  from  0  to  360  degrees  with  a  tilt 
adjustment  of  ±  10  degrees  in  azimuth  and  an  accuracy  of  0.1  degree.  The  coarse  adjustment 
in  azimuth  is  simply  achieved  by  rotating  the  whole  assembly  mounted  on  a  turn  table.  After 
reflection  on  the  scene  mirror,  the  beam  is  then  successively  focussed  by  the  Cassegrain 
telescope  and  reflected  by  an  off-axis  parabolic  mirror  (PM)  to  produce  a  collimated  beam  of 
appropriate  diameter  at  the  input-port  of  the  interferometer.  The  two  output  modules  are 
identical,  except  that  channel  1  includes  an  MCT  detector  (HgCdTe)  optimized  for  the 
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INTERFEROMETER 


FIGURE  1  -  Schematic  representation  of  the  Fourier  spectrometer  system  referred  to  as  the 
Double  input  Beam  Interferometer  Sounder  (DBIS) 
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5-20  Jim  spectral  region,  while  channel  2  includes  an  InSb  detector  optimized  for  the  2-5  pm 
region.  Note  that  the  InSb  module  is  not  explicitly  shown  in  Fig.  1.  These  modules  contain 
parabolic  and  condensing  mirrors  that  focus  the  beam  coming  from  the  interferometer  onto  a 
detector  (MCT  or  InSb)  of  1  mm  diameter.  An  aperture  wheel  (AW)  mounted  with  stops  of 
different  diameters  allows  for  the  adjustment  of  the  field  of  view  (FOV)  of  the  instrument  (5 
mrad  or  smaller).  A  camera  boresighted  to  the  telescope  can  be  used  to  aim  and  visualize  the 
measured  scene. 

The  spectrometer  is  a  standard  BOMEM  MB100  field  interferometer  that  uses  comer 
reflectors  rather  than  the  usual  flat  plate  mirrors  normally  utilized  in  the  Michelson 
configuration.  With  such  an  interferometer,  two  beams  coming  from  different  scenes  and 
entering  different  input  ports  (one  or  two)  can  be  optically  combined  and  spectrally 
subtracted  in  real  time.  Basically,  this  capability  arises  from  the  fact  that  the  two 
interferograms  associated  to  each  input  port  are  out  of  phase.  However,  in  this  study, 
measurements  were  taken  with  the  single  telescope  configuration,  where  in  this  case  the 
second  telescope  is  replaced  by  a  stable  cold  blackbody  (liquid  nitrogen;  Input  Optics-2). 

The  radiometric  calibration  of  DBIS  spectra  is  based  on  the  two-temperature  method. 
Two  temperatures  are  required  because  there  are  two  unknowns:  the  gain  and  the  offset  (self 
emission)  of  the  instrument.  To  this  end,  two  large  flat  plate  blackbodies  (12”xl2")  kept  at 
constant  temperatures  have  been  used:  one  referred  to  as  the  hot  blackbody  is  stabilized  at  a 
temperature  of  40°  C,  while  the  other  referred  to  as  the  ambient  blackbody  is  stabilized  at  a 
temperature  near  the  ambient  temperature  of  operation  (~  10  to  25°  C).  The  accuracy  of 
both  blackbodies  is  of  the  order  of  0.1°  C.  To  obtain  a  calibrated  emission  spectrum,  three 
measurements  are  required:  the  raw  spectrum  corresponding  to  the  emission  of  the 
atmosphere  and  two  raw  spectra  from  the  reference  blackbodies.  Each  raw  spectrum  is 
obtained  by  coadding  20  individual  scans  taken  at  a  resolution  of  1  cm-1.  In  these  conditions 
of  operation,  the  Root  Mean  Square  (RMS)  noise  has  been  measured  to  be  smaller  than  0.2  K 
(brightness  temperature)  almost  everywhere  in  the  usable  portion  of  the  3-20  pm  region. 
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3.0  ATMOSPHERIC  EMISSION  SPECTRA:  RESULTS 


Examples  of  downwelling  emission  spectra  (zenith  angle  =  0  degree)  expressed  in 
brightness  temperature  appear  in  Figs.  2  and  3.  The  brightness  temperature  Tg  is  simply 
obtained  by  inverting  the  Planck  radiance  equation 


Tb(v)  =  — 
In 


1.439  v 


1.191x10-12v3  ,  , 
R(v) 


where  numerical  constants  are  defined  to  match  the  radiance  spectrum  R(v)  (watts  /  cm2-ster- 
cm'1)  with  the  temperature  unit  (K)  assuming  a  wavenumber  v  in  cm-1.  The  upper  curve  of 
Fig.  2  corresponds  to  a  typical  spectral  radiance  measured  in  clear  sky  conditions.  For 
comparison  purposes,  the  bottom  curve  represents  the  corresponding  calculation  performed 
with  a  recent  version  of  the  FASCOD3  model  (Ref.  7).  Radiosonde  profiles  of  pressure, 
temperature  and  humidity  have  been  used  as  input  for  calculations  together  with  a  clear 
midlatitude  summer,  atmosphere  for  the  other  species.  Typically,  clear  sky  calculations 
adequately  reproduce  the  measurements.  Spikes  appearing  in  the  1500-1800  cm-1  region 
and  near  660  cm-1  in  the  measured  spectrum  of  Fig.  2  are  essentially  artefacts  due  to 
limitations  in  the  calibration  procedure  when  applied  to  very  strong  absorption  bands. 

The  general  shape  of  clear  sky  spectra  is  dominated  by  the  temperature  and  the 
absorption  properties  of  the  constituents  along  the  optical  path.  In  spectral  regions  where 
absorption  is  strong,  the  atmosphere  radiates  its  energy  like  a  blackbody  at  a  temperature 
roughly  equal  to  that  near  the  sounder.  This  is  observed  in  the  two  strong  absorption  bands 
of  CO2  near  15  pm  (660  cm-1)  and  4.3  pm  (2325  cm-1),  and  also  in  the  1400-1800  cm-1 
region  (6.3  pm)  where  the  absorption  is  mainly  due  to  water  vapor.  If  we  exclude  the 
important  emission  band  of  the  O3  near  9.6  pm  (1040  cm-1),  the  atmospheric  window 
region  from  800  to  1200  cm-1  (8-12  pm)  is  dominated  by  the  water  vapor  absorption  lines 
(fine  spectral  structure)  and  continuum  (smooth  spectral  shape). 

Figure  3  compares  radiances  associated  to  clear  and  cloudy  skies.  For  instance,  curve 
2  represents  typical  emission  spectra  for  a  sky  having  a  thin  cloud  of  moderate  transparency 
at  an  altitude  of  approximately  1.4  km.  In  this  case,  the  radiance  recorded  in  the  8  to  12  pm 
region  (800-1200  cm-1)  js  largely  dominated  by  the  cloud  emission,  while  in  the  3-4  pm 
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FIGURE  2  -  Typical  example  of  spectra  measured  with  the  DBIS  and  calculated  with 
FASCOD3  for  a  clear  day  corresponding  to  midlatitude  summer  conditions 
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FIGURE  3  -  Typical  ground-based  measurements  of  the  downwelling  spectral  radiance  for 
three  atmospheric  conditions:  clear  sky  (curve  1),  thin  transparent  cloud  (curve 
2),  and  opaque  cloud  (curve  3) 
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region  the  dominant  radiance  contribution  is  due  to  solar  scattering  on  cloud  particles. 
Curve  3  represents  approximately  the  same  situation  as  in  curve  2,  except  that  the  cloud  layer 
is  opaque.  In  this  case  the  solar  scattering  component  in  the  3-4  pm  region  disappears  due  to 
large  cloud  extinction  and  cloud  emission  becomes  the  dominant  contribution. 

4.0  METHOD  OF  CLQ.UD-EYALUATLQN 

Ground-based  observations  of  cloud  emissions  in  the  IR  has  prompted  the 
development  of  a  simple  spectral  method  for  estimating  the  transmittance  (at  9.3  pm)  and  the 
apparent  temperature  of  a  cloud  layer.  Figure  4  introduces  the  method  where  three  radiance 
spectra  recorded  in  the  8  to  12  pm  window  are  compared  for  different  sky  conditions:  clear, 
thin  cloud  and  opaque  cloud.  As  seen,  the  cloud  amount  tends  to  decrease  the  contrast  of  the 
ozone  emission  band  near  9.6  pm  and  the  cloud  temperature  tends  to  increase  the  overall 
emission  in  the  window.  To  take  advantage  of  this  behavior  and  retain  the  benefit  of  a  simple 
estimation  procedure,  two  narrow  spectral  channels  have  been  selected  (vertical  bars  in  Fig. 
4).  The  choice  of  this  pair  is  not  arbitrary.  Channel  1  corresponds  to  a  wavelength  of  9.48 
pm  (1054-1055  cm-1):  It  is  used  to  probe  the  ozone  signal  and  is  free  of  absorption  lines. 
Channel  2  corresponds  to  a  wavelength  of  9.15  p.m  (1093-1094  cm-1):  It  depends  essentially 
on  the  water  vapor  continuum  and  is  also  free  of  absorption  lines.  Because  the  mixing  ratio 
of  the  ozone  peaks  at  an  altitude  of  27  km,  an  important  part  of  the  ozone  emission  (channel 
1)  occurs  above  cloud  layers.  Consequently,  the  signal  received  on  channel  1  is  mainly  due 
to  the  ozone  emission  attenuated  by  the  cloud  layer.  Channel  2  may  be  seen  as  a  probe  for 
the  apparent  temperature  .  Also  shown  in  Fig.  4  are  the  absorption  coefficients  of  liquid 
water  and  ice  crystal.  These  indicate  that  water  clouds  and  ice  clouds  should  have 
approximately  the  same  optical  properties  near  the  9.3  pm  spectral  region.  In  the  11-12  pm 
region,  distinct  ice  and  water  coefficients  would  impact  differently  the  resulting  radiance. 

In  Fig.  5,  the  bulk  reflectance  of  liquid  water  and  ice  crystal  are  compared.  The  low 
reflectance  values  (~  1  %)  in  the  region  of  the  two  probing  channels  indicate  that  the 
scattering  process  inside  the  cloud  is  dominated  by  absorption  and  transmission  mechanisms. 
This  constitutes  an  advantage  in  the  proposed  method  because  it  is  assumed  that  earth  IR 
emission  reflected  down  by  the  cloud  can  be  neglected,  which  greatly  simplifies  the  equations 
of  the  method  (as  seen  below). 
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FIGURE  4  -  Spectral  radiances  for  three  different  sky  conditions  (upper  curves,  left  scale). 

Absorption  coefficients  of  liquid  water  and  ice  cristal  (bottom  curves,  right 
scale).  Vertical  bars  correspond  to  two  narrow  channels  selected  for  cloud 
temperature  and  transmittance  evaluations 


FIGURE  5  -  Bulk  reflectance  of  liquid  water  and  ice  cristal.  The  two  vertical  bars 
correspond  to  the  narrow  channels  selected  for  cloud  temperature  and 
transmittance  evaluation 
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In  the  following,  the  equations  of  the  method  are  developed  in  relation  with  the 
atmospheric  layering  depicted  in  Fig.  6.  For  a  plane-parallel  clear  atmosphere  (Fig.  6a)  in 
local  thermodynamic  equilibrium  with  no  scattering,  the  downwelling  spectral  radiance  R(v) 
may  be  expressed  in  terms  of  three  components: 

R(v)  =  ul+i(v)  tl+i  +  B(v)  [l-xm(v)]  xL(v)  +  Dl(v)  [i] 

where  L  is  an  arbitrary  level,  ul+i(v)  corresponds  to  the  downwelling  radiance  measured  at 
level  L+l,  xm(v)  is  the  molecular  transmittance  of  the  layer  L,  xL(v)  is  the  transmittance  from 
level  L  to  the  ground,  B(v)  is  the  Planck  radiance  associated  to  layer  L  and  Dl(v) 
corresponds  to  the  radiance  of  the  atmosphere  below  level  L.  Similarly,  the  downwelling 
radiance  R'(v)  of  an  atmosphere  containing  a  cloud  layer  at  level  L  (Fig.  6b)  is  given  by 

R'(v)  =  UL+l(v)  XL+l(v)  Tc  +  B(v)  [l-XcTm(V)j  XL(V)  +  DL(V)  [2] 

where  xc  is  defined  as  the  transmittance  of  the  cloud  layer,  [l-xcxm(v)]  B(v)  is  the  associated 
gray  body  radiance  of  the  cloud.  It  is  important  to  note  that  the  Planck  radiance  terms  B(v) 
in  eqs.  1  and  2  are  assumed  to  be  identical.  Also  note  that  in  eq.  2,  the  cloud  reflected 
radiation  coming  from  the  earth  emission  and  the  surrounding  atmosphere  is  neglected.  This 
represents  a  good  approximation  because  of  the  low  reflectance  of  liquid  water  and  ice  crystal 
in  the  spectral  region  near  9.3  pm.  Inserting  the  clear  sky  radiance  eq.  1  into  eq.  2  for  the 
two  probing  channels,  va  (1054-1055  cm-1)  and  Vb  (1093-1094  cm'l)  and  assuming  that  xc 
is  the  same  for  the  two  channels  yields  the  following  system  of  equations: 

R'(Va)  -  R(Va)  =  (1-XC)  [  B(va)  -  UL+l(Va)  ]  XL+l(Va)  [3] 

R’(vb)  -  R(vb)  =  (l-tc)  [  B(vb)  -  uL+i(vb)  ]  xL+i(vb) .  [4] 

These  equations  were  obtained  using  the  expression  for  the  transmittance  xL+1  =  xm  x^. 

The  two  unknowns  to  be  solved  are  xcand  B.  Other  quantities  are  properly  estimated 
based  on  clear  sky  radiance  measurements  or  calculations,  associated  radiosonde  or  inverted 
profiles  and  radiative  transfer  calculations  to  evaluate  ul+1  and  xL+1.  For  a  given 
temperature,  the  Planck  radiances  B(va)  and  B(vb)  evaluated  at  two  arbitrary  wavenumbers 
are  related  to  each  other  by  the  general  expression: 
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FIGURE  6  -  Atmospheric  layering  related  to  the  method  for  cloud  temperature  and 
transmittance  evaluation:  a-  clear  atmosphere;  b-  atmosphere  with  a  cloud  layer 
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B(vb)  = 


1  + 


1+c 


v! 

B(Va) 


[5] 


where  the  constant  c  =  1.191x1  O'12.  In  the  case  of  the  two  probing  channels  selected  for  the 
present  method,  va  (1054-1055  cm-1)  and  Vb  (1093-1094  cm-1),  eq.  5  reduces  to 

B(vb)  =  1.42198  [B(va)]  t1-0369*^  [6] 


and  eqs.  3  and  4  can  be  combined  to  eliminate  the  unknown  transmittance  xc,  which  gives  a 
simple  expression  for  B(va)  =  Ba 

Ci(Ba)  1-03698  +  C2(Ba)  +C3  =  0,  [7] 


where 


Ci  =  1.42198  mxL+1(vb)) 

C2  =  -  tL+i(va)> 

c3  =  uL+i(va)XL+i(va)  -  m  uL+i(vb)xL+1(vbX 

and 


Ba  is  found  by  solving  eq.  7  numerically  using  the  Newton  method.  In  this  case,  the  initial 
solution  B°a  is  found  by  assuming  that  Ba  =  Bb  in  eqs.  3  and  4.  The  solution  (one  iteration) 
is  given  by 


Ba  =  B«  - 


F(Bg) 

F(B2), 


[8] 


where 

f(  b2)  =  Ci  (  b2)  l03698  +  C2  ( B2)  +  c3? 
F’(  B2)  =  1.03698Ci  (  B2)  °’°3698  +  C2> 


and 
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1.42198 


+  C2 


Knowing  Ba ,  it  is  easy  to  solve  analytically  for  xc,  which  can  be  obtained  from 


_  j  |  R'a  -  R'b  -  Ra  +  Rb 


C4  -  C5 


[9] 


with 

C4  =  [  Ba  -  uL+i(va)  ]  XL+l(va) , 
and 

C5  =  [  Bb  -  UL+l(Vb)  ]  XL+l(Vb)  . 

Equations  7  and  9  summarize  the  result  of  the  cloud  evaluation  method.  The  apparent 
temperature  of  the  cloud  layer  is  simply  found  by  inverting  the  Planck  radiance  solved  at  eq. 
7.  It  depends  on  the  measured  radiances  of  the  cloudy  skies  R'a  and  R'b  ,  and  the  calculated 
(or  measured)  radiances  of  the  associated  clear  skies  Ra  and  Rb-  The  IR  transmittance  of  the 
cloud  at  9.3  |xm  is  estimated  by  eq.  9  and  depends  on  the  Planck  radiances  of  the  clouds  Ba 
and  Bb  (found  by  eq.  7),  the  measured  radiances  of  the  cloudy  skies  R'a  and  R’b  and  the 
calculated  (or  measured)  radiances  of  the  associated  clear  skies  Ra  and  Rb-  The  constants  Ci, 
C2,  C3,  C4  and  C5  are  computed  with  MODTRAN  based  on  radiosonde  or  inverted  profiles 
obtained  from  the  current  or  a  representative  clear  sky  condition. 

5.0  RESULTS  AND  ANALYSES. 

Verifications  of  the  cloud  evaluation  method  have  been  performed  on  a  case  (6 
December  1995)  in  which  coincident  radiance,  radiosonde  and  lidar  data' were  recorded.  The 
radiosonde  launch  site  was  located  at  Canadian  Forces  Base  Valcartier  (3  km  from  DREV), 
while  the  site  for  spectral  and  lidar  measurements  was  located  at  DREV.  Typically, 
radiosondes  were  launched  once  a  day  at  09:30  a.m.  local  time.  For  the  6  December  case  the 
temperature  at  ground  level  was  271  K  and  the  sky  was  totally  covered  by  a  cloud  layer  of 
variable  opacity.  Inspection  of  the  lidar  returns  indicates  a  cloud  base  at  an  altitude  of  1.4 
km  and  a  cloud  thickness  of  approximately  0.4  km.  At  this  altitude  range  (1.4-1. 8  km)  the 
radiosonde  data  exhibit  a  very  weak  temperature  inversion,  averaging  around  260  K.  This 
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corresponds  to  the  apparent  temperature  of  the  cloud  since  the  atmosphere  is  assumed  to  be 
in  thermal  equilibrium. 

5.1  Verification  of  the  Method 

A  first  verification  of  the  method  was  achieved  on  simulations.  They  were  calculated 
with  a  version  of  the  radiative  transfer  model  MODTRAN2.  For  that,  meteorological 
parameters  for  model  calculations  correspond  to  the  radiosonde  temperature  and  humidity 
profiles  (6  December  1995),  while  other  needed  profiles  for  the  dry  atmosphere  correspond 
to  midlatitude  winter  conditions.  Based  on  these  profiles,  MODTRAN2  was  run  in  two 
consecutive  steps  to  generate  all  the  required  input  to  verify  the  method  (eqs.  7  and  9).  At 
the  first  step,  MODTRAN2  was  used  to  compute  the  clear  sky  transmittance  xL+1  from  level 
L+l  to  the  ground  and  to  compute  the  downwelling  clear  sky  radiances  uL+1  and  R  received 
at  level  L+l  and  at  the  ground,  respectively.  At  the  second  step,  MODTRAN2  was  used  to 
simulate  the  downwelling  spectral  radiance  of  the  cloudy  atmosphere  R’.  Note  that  for  these 
simulations  the  multiscattering  option  was  not  utilized.  The  cirrus  cloud  model  included  in 
MODTRAN2  was  selected  because  it  offers  enough  flexibility  to  select  the  cloud  altitude  and 
thickness.  To  be  consistent  with  the  measurements,  these  parameters  were  set  at  1.4  km  and 
0.4  km,  respectively.  In  addition  to  altitude  and  thickness,  the  cirrus  cloud  option  allows  for 
the  adjustment  of  the  cloud  IR  transmittance  through  an  absorption  coefficient.  With  this 
capability,  several  simulated  radiances  have  been  generated  for  different  cloud  transmittances. 
These  simulations  were  then  used  to  test  eqs.  7  and  9. 

The  results  of  the  simulation  tests  are  gathered  in  Fig.  7.  The  reference  curve  is 
defined  as  the  "truth",  which  correspond  to  the  temperature-transmittance  input  used  in 
MODTRAN2  simulations.  For  comparison,  the  apparent  temperature  and  cloud  transmittance 
(9.3  pm)  derived  from  the  evaluation  method  are  also  shown  in  Fig.  7  (simulation  points).  As 
seen,  the  overall  agreement  between  the  two  results  is  very  good.  At  high  transmittance  level 
the  cloud  evaluation  method  slightly  differs  from  the  reference  curve  mainly  because  of 
certain  limitations  in  the  numerical  accuracy  of  MODTRAN2  calculations.  At  worst,  this 
corresponds  to  temperature  errors  smaller  than  1  degree  for  transmittances  higher  than  80  %. 
Of  less  importance  is  an  error  introduced  by  the  radiance  model  itself  such  as  formulated  in 
eq.  1.  This  equation  is  analytically  exact  only  within  the  limits  of  a  monochromatic  radiation 
(infinite  resolution).  However,  through  a  series  of  calculations,  we  have  verified  that  eq.  1  can 
also  be  used  to  compute  radiances  at  a  moderate  resolution  of  1  cm'1  where  the  resulting 
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FIGURE  7  -  Results  of  the  cloud  temperature-transmittance  estimation  method  for  different 
simulated  and  measured  spectra  corresponding  to  the  6  December  95 
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FIGURE  8  -  Sequence  of  spectral  radiance  measurements  in  cloudy  sky  conditions  recorded 
the  6  December  95  between  9:19  a.m.  and  10:05  a.m. 
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numerical  error  appears  to  be  negligible. 

Secondly,  the  cloud  evaluation  method  has  been  tested  on  a  sequence  of  experimental 
spectra  reported  in  Fig.  8.  These  spectral  radiance  measurements  correspond  to  a  decreasing 
cloudiness  recorded  on  December  6th,  1995  between  9:19  a.m.  and  10:05  a.m.  The 
corresponding  temperature  and  transmittance  results  appear  in  Fig.  7.  In  this  case,  the 
evaluation  method  gives  the  right  order  of  magnitude,  but  there  are  differences  of  up  to  6  K. 
The  temporal  variability  of  the  cloud  parameters  with  respect  to  the  lidar  sounding 
(height=1.4  km  and  thickness=0.4  km)  is  partly  responsible  for  this  discrepancy:  For  the 
time-coincident  lidar  measurement  (round  dot  in  Fig.  7),  the  temperature  difference  is  smaller 
than  2  K.  A  remaining  source  of  error  may  also  be  due  to  the  fact  that  the  actual  line  of  sight 
of  both  sounding  systems,  the  spectrometer  and  the  lidar,  was  not  exactly  coincident.  In  this 
case  the  spatial  variability  of  the  cloud  parameters  may  be  the  limiting  factor.  Unfortunately, 
with  the  actual  setup,  it  was  not  possible  to  aim  at  exactly  the  same  cloud  area  with  the  two 
instruments. 

Finally,  as  stated  earlier,  a  minor  source  of  errors  comes  from  the  fact  that  the 
radiance  measurements  contain  an  additional  contribution  from  the  multiply  scattered 
radiation  in  the  cloud,  while  the  evaluation  method  ignores  this  contribution.  For  the  time 
being,  there  is  no  reliable  radiative  transfer  model  available  to  properly  evaluate  the  multi¬ 
scattering  effect  in  cloudy  atmosphere:  Recent  versions  of  models  such  as  MODTRAN  and 
FASCODE  offers  multiscattering  algorithm,  but  they  have  not  been  fully  validated  yet. 
However,  simple  arguments  suggest  that  this  multiscattering  contribution  would  be  negligible 
in  the  9.3  pm  region  because  the  reflectance  of  cloud  constituents  (liquid  water  and  ice 
crystal)  is  very  small  (approximately  1%,  see  Fig.  5),  which  tends  to  favor  absorption  rather 
than  scattering  in  the  extinction  process. 

5.2  Extension  of  the  Method 

The  multiscattering  effects  and  the  parasite  contribution  from  the  earth  radiation 
reflected  from  the  cloud  such  as  discussed  in  Section  4,  remain  two  issues  that  should  be 
addressed  in  further  model  developments.  With  a  proper  handling  of  these  two  radiance 
contributions  it  would  be  possible  to  extend  the  current  model  to  the  evaluation  of  cloud 
transmittance  in  the  10  to  12  pm  region  (800-1000  cm-1).  The  proposed  procedure  of  cloud 
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evaluation  in  the  context  of  an  extended  method  can  be  summarized  in  the  six  following 
steps: 

1-  Evaluation  of  the  cloud  temperature  and  transmittance  by  eqs.  7  and  9, 

2-  Inversion  of  the  lidar  data  to  extract  the  mean  radius  of  the  cloud  particles  and  the 
extinction  coefficient, 

3-  Computation  of  both  the  multiscattering  and  the  cloud-reflected  radiance  contributions 
(8-12  pm)  from  the  inverted  cloud  parameters; 

4-  Correction  of  the  measured  spectral  radiance  to  remove  these  two  contributions  (cloud- 
reflected  and  multiscattering); 

5-  Second  application  of  eqs.  7  and  9  on  the  corrected  radiance  spectrum,  but  with  an 
extension  to  other  micro-windows  found  in  the  IR.  This  would  allow  for  the 
determination  of  the  cloud  transmittance  in  the  800-1000  cm'l  region  (8-10  pm).  In  this 
spectral  domain,  distinct  optical  properties  of  ice  crystal  and  liquid  water  make 
transmittance  sensitive  to  the  nature  of  the  cloud  particles.  In  this  case,  the  water  /  ice 
cloud  ratio  can  be  obtained  from  the  shape  of  the  spectral  transmittance; 

6-  Computation  of  more  accurate  cloud-reflected  and  multiscattering  radiance  contributions 
using  the  results  of  step  5,  generation  of  an  updated  correction  for  the  measured  radiance, 
and  again  application  of  eqs.  7  and  9  to  improve  the  cloud  temperature-transmittance 
determination. 

The  benefit  of  this  iterative  method  of  cloud  evaluation  based  on  spectral  and  lidar  data  is  to 
provide  all  the  microphysical  information  needed  to  model  the  radiative  transfer  throughout 
the  visible  up  to  the  far  infrared. 

A  fundamental  application  of  this  extended  method  is  to  help  solving  the  often 
reported  ambiguity  concerning  daytime  radiance  in  the  1  to  4  jam  region.  In  this  case,  the 
strong  solar  scattering  on  the  cloud  particles  is  poorly  modelled  by  most  of  the  available 
radiative  transfer  models.  Factors  of  two  between  simulated  and  measured  radiances  are  not 
uncommon  in  this  spectral  region.  This  is  one  of  the  main  concerns  of  today’s  modelers.  It 
is  believed  that  the  method  of  cloud  evaluation  based  on  the  spectral  and  lidar  data  proposed 
in  this  report  might  constitute  a  powerful  tool  to  support  the  development  and  validation  of 
radiative  transfer  models  for  cloudy  atmospheres. 
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5.3  Impa.ct.-Qf, Background. Profiles 

An  important  aspect  that  needs  to  be  addressed  is  the  impact  of  background  profiles 
used  in  the  method  for  the  calculations  of  clear  sky  radiance  and  transmittance  data  (R,  ul+i 
and  TL+l)-  Figure  9  shows  radiosonde  recordings  of  temperature  and  humidity  profiles  for 
December  5th  (clear  day)  and  6th  (cloudy  day),  1995.  It  is  convenient  to  remember  that  the 
temperature- transmittance  results  shown  in  Fig.  7  were  derived  based  on  the  exact  profiles  of 
the  day  of  measurements  (6  Dec.  95).  The  same  format  of  results,  but  computed  with 
approximate  profiles  (typical)  rather  than  the  exact  ones  (actual),  appears  in  Fig.  10.  The 
approximate  profiles  are  simply  those  recorded  the  day  before  (5  Dec.  95)  in  clear  sky 
conditions.  As  seen,  the  method  is  weakly  sensitive  to  the  input  profiles.  Even  with  errors  of 
10  K  in  the  first  kilometer  for  the  temperature  profile  together  with  a  level  of  humidity  two 
times  lower  than  the  actual  one,  the  induced  errors  in  the  cloud  temperatures  and 
transmittances  are  approximately  2  K  and  4%,  respectively. 

This  low  sensitivity  to  the  input  profiles  has  the  practical  advantage  that  the  cloud 
evaluation  method  can  be  utilized  without  supporting  radiosondes.  The  previous  example 
corresponds  to  an  extreme  scenario  in  which  no  additional  effort  was  made  to  improve  the 
profile  estimates.  In  real  situations,  the  ground  level  temperature  and  humidity  are  easily 
available.  Thus,  the  application  of  a  simple  technique  for  scaling  the  rough  profiles  with 
ground  level  meteorological  data  would  result  in  improved  profile  estimates.  The 
temperature  and  humidity  profiles  can  also  be  improved  based  on  the  radiance  spectrum 
itself.  For  instance,  in  Ref.  6,  a  simple  algorithm  has  been  proposed  to  estimate  the 
temperature  profiles  in  the  first  500  meters  (above  ground)  using  the  spectral  radiance  in  the 
CO2  absorption  band  near  15  pm.  Another  practical  example  is  the  case  of  broken  clouds. 
In  this  case,  the  downwelling  spectral  radiance  from  a  clear  portion  of  the  sky  can  be  inverted 
to  retrieve  the  full  temperature  and  humidity  profiles  as  proposed  in  (Refs.  2  and  7).  It  is 
believed  that,  in  operational  scenarios  and  without  a  radiosonde,  the  systematic  use  of  a  simple 
technique  for  obtaining  proper  profile  estimates  can  improve  the  cloud  evaluation  method  to 
a  level  of  approximately  1  K  error  in  temperature  and  2%  in  transmittance. 

To  further  illustrate  the  potential  of  the  method,  Figs.  11  and  12  show  another 
example  in  which  the  method  has  been  applied  to  a  sequence  of  radiance  spectra  recorded  in 
the  cloudy  conditions  of  December  4th,  1995  (Fig.  11).  There  was  no  radiosonde  available 
for  this  day.  However,  the  ground  temperature  and  humidity  were  very  similar  to  the 
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FIGURE  9  -  Profiles  of  temperature  and  absolute  humidity  recorded  (9:30  a.m.)  by  the 
radiosonde  system  for  December  5th  (clear  sky)  and  6th  (cloudy  sky),  1995 
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FIGURE  10  -  Impact  of  the  background  profiles:  The  cloud  temperature-transmittance 
estimation  method  is  applied  using  typical  profiles  (5  December  1995)  rather 
than  the  actual  ones  (6  December  1995)  . 
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Sequence  of  spectral  radiance  measurements  in  cloudy  sky  conditions 
recorded  on  December  4th,  1995  between  9:56  a.m.  and  10:18  a.m. 
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FIGURE  12  -  Cloud  temperature-transmittance  estimations  for  a  sequence  of  spectra 
recorded  on  December  4th,  1995 


P500380.PDF  [Page:  30  of  44] 


UNCLASSIFIED 

20 


ones  recorded  the  following  day  (5  Dec.  95).  For  this  reason,  the  radiosonde  temperature 
and  humidity  profiles  available  for  Dec.  5^  were  utilized  as  a  good  estimate  to  process  the 
cloudy  sky  spectra  of  Dec.  4th.  The  results  of  the  cloud  evaluation  are  depicted  in  Fig.  12. 
As  seen,  for  each  spectrum  of  the  sequence  (9:56  a.m.  to  10:18  a.m.),  the  cloud  temperature 
is  very  stable  and  close  to  259  K  on  the  average,  while  transmittance  varies  by  approximately 
a  factor  of  two  from  25%  to  55%.  Unfortunately,  because  of  no  radiosonde  was  available  for 
this  day,  it  was  not  possible  to  directly  validate  the  cloud  temperature  evaluations. 

However,  in  a  preliminary  study,  the  current  multiscattering  inversion  algorithm  of 
Bissonnette  (Ref.  1)  has  been  applied  to  the  lidar  data  recorded  concurrently  with  the  spectral 
data.  The  result  of  the  lidar  inversion  suggests  that  the  average  radius  of  the  cloud  particles  is 
roughly  equal  to  2  pm,  which  leads  to  a  cloud  transmittance  of  approximately  2  %  at  the  lidar 
wavelength  (1.06  pm),  This  value  of  transmittance  is  in  good  agreement  with  the  cloud 
opacity  noted  by  visual  observations.  It  also  seems  to  contrast  the  high  transmittance  level 
found  in  the  IR  by  the  cloud  evaluation  method  (Fig.  12).  However,  using  simple  arguments 
borrowed  from  the  diffraction  theory  applied  to  spherical  particles,  it  appears  that  the  high 
values  of  cloud  transmittance  derived  from  spectral  data  at  9.3  pm  (25-55%)  are  very 
consistent  with  cloud  particle  radius  of  2  pm  i.  e.  cloud  transmittance  results  are  in  good 
agreement  with  what  is  expected  from  lidar  data. 

5.4  Impact  of  Ozone  Variability 

The  last  major  source  of  possible  limitations  in  the  implementation  of  the  cloud 
evaluation  method  is  the  variability  of  ozone  concentration.  For  now,  the  method  is  based  on 
the  spectral  contrast  of  the  ozone  emission  near  9.  6  pm.  Thus,  to  be  accurate,  it  requires 
good  knowledge  of  the  ozone  concentration.  It  is  well  known  that  the  ozone  profiles  in  the 
atmosphere  are  quite  stable  from  time  to  time  but  they  vary  with  the  location  and  season. 
Usually,  the  ozone  mixing  ratios  are  properly  accounted  for  by  the  six  standard  profiles: 
tropical,  midlatitude  summer,  midlatitude  winter,  subarctic  summer,  subarctic  winter  and  U.  S. 
standard.  These  are  the  profiles  included  in  MODTRAN  that  should  be  used  in  relation  to  the 
proposed  cloud  evaluation  method.  However,  the  industrial  activities  in  and  near  the  urban 
areas  favor  both  the  creation  and  destruction  of  ozone  near  the  surface:  Punctual 
measurements  performed  in  polluted  areas  usually  show  variations  of  ozone  concentration  as 
important  as  5  to  60  ppm  over  a  period  of  a  few  hours  with  a  repetitive  cycle  of 
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FIGURE  13 
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Impact  of  ozone  variability  near  the  surface  :  a-  the  O3 -concentration  is 
reduced  to  a  third  of  the  actual  value  in  the  first  500  m  above  the  ground,  while 
for  b-  it  is  increased  to  three  times  the  actual  value  (<  500  m).  For  these 
calculations,  the  O3  profile  is  kept  equal  to  the  default  value  for  this  location 
and  season  (midlatitude  winter). 
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approximately  24  hours.  These  measurements  are  usually  taken  very  close  to  the  ground 
(approximately  3  m  above),  and  to  the  best  of  our  knowledge  the  characterization  of  the 
ozone  variability  in  the  first  kilometer  above  the  surface  has  not  been  addressed  yet. 

In  order  to  estimate  the  impact  of  ozone  variability  near  the  surface,  we  have 
simulated  the  effect  of  erroneous  ozone  concentrations.  Figure  13a  represents  the 
temperature-transmittance  results  obtained  when  the  O3  concentration  is  reduced  to  a  third  of 
an  assumed  typical  value  of  30  ppm  in  the  first  500  meters  above  the  ground.  In  Fig,  13b, 
the  ozone  concentration  is  three  times  the  actual  one.  Overall,  the  impact  is  much  more 
pronounced  at  high  values  of  cloud  transmittance.  For  cloud  transmittances  lower  than 
approximately  50%  the  impact  is  not  too  severe:  a  66%  error  in  the  ozone  concentration 
(near  the  surface)  yields  cloud  temperature  errors  smaller  than  2  K  and  cloud  transmittance 
errors  of  the  order  of  3%. 

Finally,  to  cover  all  possible  scenarios  of  cloud  characterization,  including  urban 
areas,  it  might  be  requested  to  provide  the  method  with  information  from  a  punctual  sensor  of 
ozone.  However,  over  the  last  three  years  many  clear  sky  measurements  of  the  ozone 
emission  band  (near  9.6  pm)  taken  at  DREV  with  the  DBIS  do  not  exhibit  such  large  ozone 
variations  as  those  recorded  by  the  punctual  sensors.  This  suggests  that  the  ozone 
fluctuations  near  the  surface  appear  in  a  layer  thickness  very  much  smaller  than  500  m,  which 
in  this  case  should  not  seriously  limit  the  accuracy  of  the  cloud  evaluation  method. 


6.0  SUMMARY  AND  CONCLUSIONS 

A  method  for  estimating  the  cloud  temperature  and  transmittance  from  the 
downwelling  spectral  radiance  has  been  proposed.  The  development  of  the  method  is 
supported  by  ground-based  observations  performed  with  a  Fourier  transform  spectrometer 
system  having  a  spectral  resolution  of  1  cm-1  and  covering  the  IR  band  from  3  to  20  pm. 

The  basic  principle  of  the  method  is  simple.  It  uses  of  the  strong  emission  band  of 
ozone  in  the  infrared.  Because  most  of  the  ozone  emission  originates  from  the  upper 
atmosphere,  it  constitutes  a  natural  source  of  IR  radiation  for  probing  clouds  in  the  lower  part 
of  the  atmosphere  using  a  ground-based  receiver.  Temperature  and  transmittance  are  derived 
from  the  cloudy  sky  radiance  measured  in  two  narrow  spectral  channels  :  one  is  in  the  ozone 
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emission  band  at  9.48  pm  and  the  second  one  is  off  the  band  at  9.15  pm.  The  choice  of 
these  channels  simplifies  the  calculation  because  they  are  free  of  water  vapor  absorption  lines. 
The  solution  for  cloud  temperature  and  transmittance  is  found  by  solving  a  system  of  two 
equations  which  requires,  as  input,  an  estimate  of  the  cloud  height  and  thickness  obtained 
from  lidar  returns,  and  an  estimate  of  atmospheric  temperature  and  humidity  profiles  to 
simulate  the  background  atmosphere  from  MODTRAN  calculations. 

The  method  has  been  tested  in  two  different  ways.  The  first  validation  was  done  to 
verify  the  impact  of  the  few  assumptions  behind  the  method.  For  that  a  series  of  simulations 
were  performed  with  the  atmospheric  model  MODTRAN2  using  the  cirrus  cloud  option. 
These  simulations  indicate  that  temperature  errors  introduced  by  the  approximations  are 
always  smaller  than  1  K  for  any  cloud  transmittance.  The  second  validation  was  done  on  a 
selected  series  of  experimental  spectra.  In  this  case,  the  evaluation  method  gives  quite 
accurate  results  with  typical  temperature  errors  smaller  than  2  K  for  spectral  measurements 
coincident  in  time  with  the  lidar  measurements.  The  main  part  of  the  remaining  error  has 
been  attributed  to  the  imperfect  boresighting  of  the  spectrometer  and  lidar  systems. 

The  study  is  completed  by  an  analysis  of  the  two  main  sources  of  limitation  that 
might  be  encountered  in  the  application  of  the  method.  The  first  limitation  is  related  to  the 
accuracy  of  background  temperature  and  humidity  profiles  required  by  the  method.  Results 
of  this  analysis  indicate  a  weak  sensitivity  to  input  profiles:  Even  with  errors  of  10  K  in  the 
first  kilometer  of  the  temperature  profile,  coupled  to  an  error  of  a  factor  of  two  in  the 
humidity  profile  the  resulting  errors  in  cloud  temperature  and  transmittance  are  usually 
smaller  than  2  K  and  4%,  respectively.  This  low  sensitivity  to  input  profiles  has  the  great 
practical  advantage  that  the  cloud  evaluation  method  can  be  utilized  without  supporting 
radiosondes.  The  second  source  of  potential  limitation  is  due  to  the  time  variability  of  ozone 
concentration.  Through  a  series  of  simulations  the  impact  of  the  ozone  variability  in  the  first 
500  m  of  the  atmosphere  above  ground  has  been  found  to  be  not  too  severe:  For  cloud 
transmittance  lower  than  50  %  an  error  of  66  %  in  ozone  concentration  near  the  surface  leads 
to  temperature  errors  smaller  than  2  K  and  cloud  transmittance  errors  of  the  order  of  3  %. 
This  level  of  accuracy  can  be  further  improved  by  adding  a  punctual  ozone  sensor  to  the 
actual  instrumentation:  spectrometer  and  lidar  systems. 

Next  activities  will  be  oriented  to  generalizing  the  cloud  evaluation  method  to  any 
slant  path  from  the  ground  or  from  an  airborne  platform.  A  particular  emphasis  will  be  put 
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on  evaluating  the  impact  of  multiple  scattering  on  cloud  temperature  and  transmittance 
determinations.  An  alternate  method  which  uses  two  pure  water  vapor  continuum 
microwindows  to  get  ride  of  the  ozone  variability  constraint  will  be  investigated. 
Theoretical  refinements  of  the  method  are  also  required  to  support  wider  spectral  channels 
potentially  usable  in  a  simplified  portable  instrument,  a  two-band  cloud  radiometer. 

This  report  provides  a  basis  to  further  develop  the  cloud  evaluation  method  for 
applications  in  icing  hazard  detection.  The  ultimate  goal  is  to  demonstrate  that  combined 
lidar  and  spectral  measurements  can  be  used  for  determining  the  temperature,  the 
concentration  and  the  composition  (ice  crystals  or  supercooled  water  droplets)  of  clouds. 
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